Comparison of literature mining tools for variant classification: Through the lens of 50 RYR1 variants.

PURPOSE: The American College of Medical Genetics and Genomics and the Association for Molecular Pathology have outlined a schema that allows for systematic classification of variant pathogenicity. Although gnomAD is generally accepted as a reliable source of population frequency data and ClinGen has provided guidance on the utility of specific bioinformatic predictors, there is no consensus source for identifying publications relevant to a variant. Multiple tools are available to aid in the identification of relevant variant literature, including manually curated databases and literature search engines. We set out to determine the utility of 4 literature mining tools used for ascertainment to inform the discussion of the use of these tools. METHODS: Four literature mining tools including the Human Gene Mutation Database, Mastermind, ClinVar, and LitVar 2.0 were used to identify relevant variant literature for 50 RYR1 variants. Sensitivity and precision were determined for each tool. RESULTS: Sensitivity among the 4 tools ranged from 0.332 to 0.687. Precision ranged from 0.389 to 0.906. No single tool retrieved all relevant publications. CONCLUSION: At the current time, the use of multiple tools is necessary to completely identify the literature relevant to curate a variant.

Genome-wide identification of dominant polyadenylation hexamers for use in variant classification.

Polyadenylation is an essential process for the stabilization and export of mRNAs to the cytoplasm and the polyadenylation signal hexamer (herein referred to as hexamer) plays a key role in this process. Yet, only 14 Mendelian disorders have been associated with hexamer variants. This is likely an under-ascertainment as hexamers are not well defined and not routinely examined in molecular analysis. To facilitate the interrogation of putatively pathogenic hexamer variants, we set out to define functionally important hexamers genome-wide as a resource for research and clinical testing interrogation. We identified predominant polyA sites (herein referred to as pPAS) and putative predominant hexamers across protein coding genes (PAS usage >50% per gene). As a measure of the validity of these sites, the population constraint of 4532 predominant hexamers were measured. The predominant hexamers had fewer observed variants compared to non-predominant hexamers and trimer controls, and CADD scores for variants in these hexamers were significantly higher than controls. Exome data for 1477 individuals were interrogated for hexamer variants and transcriptome data were generated for 76 individuals with 65 variants in predominant hexamers. 3' RNA-seq data showed these variants resulted in alternate polyadenylation events (38%) and in elongated mRNA transcripts (12%). Our list of pPAS and predominant hexamers are available in the UCSC genome browser and on GitHub. We suggest this list of predominant hexamers can be used to interrogate exome and genome data. Variants in these predominant hexamers should be considered candidates for pathogenic variation in human disease, and to that end we suggest pathogenicity criteria for classifying hexamer variants.

Clinical and molecular heterogeneity of syndromic hypothalamic hamartoma.

Two children are presented who have a distinct syndrome of multiple buccolingual frenula, a stiff and short fifth finger with small nails, a hypothalamic hamartoma, mild to moderate neurological impairment, and mild endocrinological symptoms. No variant assessed to be pathogenic or likely pathogenic was detected in the GLI3 gene in either child. This syndrome appears to be distinct from the inherited Pallister-Hall syndrome associated with GLI3 variants, which is characterized by hypothalamic hamartoma, mesoaxial polydactyly, and other anomalies. In the individuals described here, manifestations outside of the central nervous system were milder and the mesoaxial polydactyly, which is common in individuals with Pallister-Hall syndrome, was absent. Instead, these children had multiple buccolingual frenula together with the unusual appearance of the fifth digit. It remains unclear whether these two individuals represent a separate nosologic entity or if they represent a milder manifestation of one of the more severe syndromes associated with a hypothalamic hamartoma.

Co-occurrence of Proteus syndrome and ventricular tachycardia cardiac arrest in a teenager.

Proteus syndrome is an extremely rare overgrowth condition caused by a somatic variant of the AKT1 gene. It can involve multiple organ systems though rarely is there symptomatic cardiac involvement. Fatty infiltration of the myocardium has been described but has not been reported to cause functional or conduction abnormalities. We present an individual with Proteus syndrome who suffered a sudden cardiac arrest.

Updated variant curation expert panel criteria and pathogenicity classifications for 251 variants for RYR1-related malignant hyperthermia susceptibility.

The ClinGen malignant hyperthermia susceptibility (MHS) variant curation expert panel specified the American College of Medical Genetics and Genomics/Association of Molecular Pathologists (ACMG/AMP) criteria for RYR1-related MHS and a pilot analysis of 84 variants was published. We have now classified an additional 251 variants for RYR1-related MHS according to current ClinGen standards and updated the criteria where necessary. Criterion PS4 was modified such that individuals with multiple RYR1 variants classified as pathogenic (P), likely pathogenic (LP), or variant of uncertain significance (VUS) were not considered as providing evidence for pathogenicity. Criteria PS1 and PM5 were revised to consider LP variants at the same amino-acid residue as providing evidence for pathogenicity at reduced strength. Finally, PM1 was revised such that if PS1 or PM5 are used PM1, if applicable, should be downgraded to supporting. Of the 251 RYR1 variants, 42 were classified as P/LP, 16 as B/LB, and 193 as VUS. The primary driver of 175 VUS classifications was insufficient evidence supporting pathogenicity, rather than evidence against pathogenicity. Functional data supporting PS3/BS3 was identified for only 13 variants. Based on the posterior probabilities of pathogenicity and variant frequencies in gnomAD, we estimated the prevalence of individuals with RYR1-related MHS pathogenic variants to be between 1/300 and 1/1075, considerably higher than current estimates. We have updated ACMG/AMP criteria for RYR1-related MHS and classified 251 variants. We suggest that prioritization of functional studies is needed to resolve the large number of VUS classifications and allow for appropriate risk assessment. RYR1-related MHS pathogenic variants are likely to be more common than currently appreciated.

Clinical exome sequencing of 1000 families with complex immune phenotypes: Toward comprehensive genomic evaluations.

BACKGROUND: Prospective genetic evaluation of patients at this referral research hospital presents clinical research challenges. OBJECTIVES: This study sought not only a single-gene explanation for participants' immune-related presentations, but viewed each participant holistically, with the potential to have multiple genetic contributions to their immune phenotype and other heritable comorbidities relevant to their presentation and health. METHODS: This study developed a program integrating exome sequencing, chromosomal microarray, phenotyping, results return with genetic counseling, and reanalysis in 1505 individuals from 1000 families with suspected or known inborn errors of immunity. RESULTS: Probands were 50.8% female, 71.5% were ≥18 years, and had diverse immune presentations. Overall, 327 of 1000 probands (32.7%) received 361 molecular diagnoses. These included 17 probands with diagnostic copy number variants, 32 probands with secondary findings, and 31 probands with multiple molecular diagnoses. Reanalysis added 22 molecular diagnoses, predominantly due to new disease-gene associations (9 of 22, 40.9%). One-quarter of the molecular diagnoses (92 of 361) did not involve immune-associated genes. Molecular diagnosis was correlated with younger age, male sex, and a higher number of organ systems involved. This program also facilitated the discovery of new gene-disease associations such as SASH3-related immunodeficiency. A review of treatment options and ClinGen actionability curations suggest that at least 251 of 361 of these molecular diagnoses (69.5%) could translate into ≥1 management option. CONCLUSIONS: This program contributes to our understanding of the diagnostic and clinical utility whole exome analysis on a large scale.

Erratum: Orofacial overgrowth with peripheral nerve enlargement and perineuriomatous pseudo-onion bulb proliferations is part of the PIK3CA-related overgrowth spectrum.

[This corrects the article DOI: 10.1016/j.xhgg.2020.100009.].

The ACMG SF v3.0 gene list increases returnable variant detection by 22% when compared with v2.0 in the ClinSeq cohort.

PURPOSE: The American College of Medical Genetics and Genomics (ACMG) recommends the return of pathogenic and likely pathogenic (P/LP) secondary findings from exome and genome sequencing. The latest version (ACMG secondary finding [SF] v3.0) includes 14 additional genes. We interrogated the ClinSeq cohort for variants in these genes to determine the additional yield in unselected individuals. METHODS: Exome data from 1473 individuals (60% White, 34% African American or Black, 6% other) were analyzed. We restricted our analyses to coding variants; +1,+2,-1, and -2 splice site variants; and the pathogenic GAA variant, NM_000152.5:c.-32-13T>G. Variants were assessed with slightly modified ACMG/Association of Molecular Pathology guidelines. RESULTS: A total of 25 P/LP variants were identified. In total, 7 individuals had P/LP variants in genes recommended for return of heterozygous variants, namely HNF1A (1), PALB2 (3), TMEM127 (1), and TTN (2). In total, 4 individuals had a homozygous variant in a gene recommended for biallelic variant return, namely HFE, NM_000410.3(HFE):c.845G>A p.Cys282Tyr. A total of 17 P/LP variants were identified in the heterozygous state in genes recommended only for biallelic variant reporting and were not returned. The frequency of returnable P/LP variants did not significantly differ by race. CONCLUSION: Using the ACMG SF v3.0, the returnable P/LP variant frequency increased in the ClinSeq cohort by 22%, from 3.4% (n = 50, ACMG SF v2.0) to 4.1% (n = 61, ACMG SF v3.0).

Somatic PIK3R1 variation as a cause of vascular malformations and overgrowth.

PURPOSE: Somatic activating variants in the PI3K-AKT pathway cause vascular malformations with and without overgrowth. We previously reported an individual with capillary and lymphatic malformation harboring a pathogenic somatic variant in PIK3R1, which encodes three PI3K complex regulatory subunits. Here, we investigate PIK3R1 in a large cohort with vascular anomalies and identify an additional 16 individuals with somatic mosaic variants in PIK3R1. METHODS: Affected tissue from individuals with vascular lesions and overgrowth recruited from a multisite collaborative network was studied. Next-generation sequencing targeting coding regions of cell-signaling and cancer-associated genes was performed followed by assessment of variant pathogenicity. RESULTS: The phenotypic and variant spectrum associated with somatic variation in PIK3R1 is reported herein. Variants occurred in the inter-SH2 or N-terminal SH2 domains of all three PIK3R1 protein products. Phenotypic features overlapped those of the PIK3CA-related overgrowth spectrum (PROS). These overlapping features included mixed vascular malformations, sandal toe gap deformity with macrodactyly, lymphatic malformations, venous ectasias, and overgrowth of soft tissue or bone. CONCLUSION: Somatic PIK3R1 variants sharing attributes with cancer-associated variants cause complex vascular malformations and overgrowth. The PIK3R1-associated phenotypic spectrum overlaps with PROS. These data extend understanding of the diverse phenotypic spectrum attributable to genetic variation in the PI3K-AKT pathway.

Low-level variant calling for non-matched samples using a position-based and nucleotide-specific approach.

BACKGROUND: The widespread use of next-generation sequencing has identified an important role for somatic mosaicism in many diseases. However, detecting low-level mosaic variants from next-generation sequencing data remains challenging. RESULTS: Here, we present a method for Position-Based Variant Identification (PBVI) that uses empirically-derived distributions of alternate nucleotides from a control dataset. We modeled this approach on 11 segmental overgrowth genes. We show that this method improves detection of single nucleotide mosaic variants of 0.01-0.05 variant allele fraction compared to other low-level variant callers. At depths of 600 × and 1200 ×, we observed > 85% and > 95% sensitivity, respectively. In a cohort of 26 individuals with somatic overgrowth disorders PBVI showed improved signal to noise, identifying pathogenic variants in 17 individuals. CONCLUSION: PBVI can facilitate identification of low-level mosaic variants thus increasing the utility of next-generation sequencing data for research and diagnostic purposes.

Clinical diagnosis of presumed SOX2 gonadosomatic mosaicism.

Purpose: To describe a family with presumed SOX2 gonadosomatic mosaicism diagnosed upon ophthalmic examination of the proband's mother.Methods: The family underwent comprehensive ophthalmic and physical examination. Variant detection was performed using trio exome analysis on peripheral leukocyte DNA from blood and saliva samples. Variant segregation analysis was performed using a custom panel NGS sequencing. An identified variant in the SOX2 gene was confirmed in the proband by Sanger sequencing.Results: We report an individual with bilateral microphthalmia, developmental delay, hearing loss, and dysmorphic features. Her mother was found to have asymptomatic forme fruste uveal coloboma affecting her anterior segment. Her father, aunt, and sisters were unaffected. Trio exome sequence analysis showed an apparent de novo heterozygous deletion in the proband, NM_003106.3:c.70_89del, NP_003097.1:p.(Asn24Argfs*65), classified as pathogenic. Testing of the other family members' peripheral blood and saliva was negative for this variant. The iris transillumination abnormalities in the proband's mother supports a gonadosomatic mosaicism scenario.Conclusions: The results from this family underscore the importance of performing detailed evaluations of the parents of apparently sporadically affected individuals with heritable ophthalmic disorders. The identification of mildly affected individuals could substantially alter recurrence risks.

Variant curation expert panel recommendations for RYR1 pathogenicity classifications in malignant hyperthermia susceptibility.

PURPOSE: As a ClinGen Expert Panel (EP) we set out to adapt the American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) pathogenicity criteria for classification of RYR1 variants as related to autosomal dominantly inherited malignant hyperthermia (MH). METHODS: We specified ACMG/AMP criteria for variant classification for RYR1 and MH. Proposed rules were piloted on 84 variants. We applied quantitative evidence calibration for several criteria using likelihood ratios based on the Bayesian framework. RESULTS: Seven ACMG/AMP criteria were adopted without changes, nine were adopted with RYR1-specific modifications, and ten were dropped. The in silico (PP3 and BP4) and hotspot criteria (PM1) were evaluated quantitatively. REVEL gave an odds ratio (OR) of 23:1 for PP3 and 14:1 for BP4 using trichotomized cutoffs of ≥0.85 (pathogenic) and ≤0.5 (benign). The PM1 hotspot criterion had an OR of 24:1. PP3 and PM1 were implemented at moderate strength. Applying the revised ACMG/AMP criteria to 44 recognized MH variants, 29 were classified as pathogenic, 13 as likely pathogenic, and 2 as variants of uncertain significance. CONCLUSION: Curation of these variants will facilitate classification of RYR1/MH genomic testing results, which is especially important for secondary findings analyses. Our approach to quantitatively calibrating criteria is generalizable to other variant curation expert panels.

Variant Classification Concordance using the ACMG-AMP Variant Interpretation Guidelines across Nine Genomic Implementation Research Studies.

Harmonization of variant pathogenicity classification across laboratories is important for advancing clinical genomics. The two CLIA-accredited Electronic Medical Record and Genomics Network sequencing centers and the six CLIA-accredited laboratories and one research laboratory performing genome or exome sequencing in the Clinical Sequencing Evidence-Generating Research Consortium collaborated to explore current sources of discordance in classification. Eight laboratories each submitted 20 classified variants in the ACMG secondary finding v.2.0 genes. After removing duplicates, each of the 158 variants was annotated and independently classified by two additional laboratories using the ACMG-AMP guidelines. Overall concordance across three laboratories was assessed and discordant variants were reviewed via teleconference and email. The submitted variant set included 28 P/LP variants, 96 VUS, and 34 LB/B variants, mostly in cancer (40%) and cardiac (27%) risk genes. Eighty-six (54%) variants reached complete five-category (i.e., P, LP, VUS, LB, B) concordance, and 17 (11%) had a discordance that could affect clinical recommendations (P/LP versus VUS/LB/B). 21% and 63% of variants submitted as P and LP, respectively, were discordant with VUS. Of the 54 originally discordant variants that underwent further review, 32 reached agreement, for a post-review concordance rate of 84% (118/140 variants). This project provides an updated estimate of variant concordance, identifies considerations for LP classified variants, and highlights ongoing sources of discordance. Continued and increased sharing of variant classifications and evidence across laboratories, and the ongoing work of ClinGen to provide general as well as gene- and disease-specific guidance, will lead to continued increases in concordance.

Factors associated with provision of smoking cessation support to pregnant women - a cross-sectional survey of midwives in New South Wales, Australia.

BACKGROUND: Smoking is the most important preventable cause of adverse pregnancy outcomes, but provision of smoking cessation support (SCS) to pregnant women is poor. We examined the association between midwives' implementation of SCS (5As - Ask, Advise, Assess, Assist, Arrange follow-up) and reported barriers/enablers to implementation. METHODS: On-line anonymous survey of midwives providing antenatal care in New South Wales (NSW), Australia, assessing provision of the 5As and barriers/enablers to their implementation, using the Theoretical Domains Framework (TDF). Factor analyses identified constructs underlying the 5As; and barriers/enablers. Multivariate general linear models examined relationships between the barrier/enabler factors and the 5As factors. RESULTS: Of 750 midwives invited, 150 (20%) participated. Respondents more commonly reported Asking and Assessing than Advising, Assisting, or Arranging follow-up (e.g. 77% always Ask smoking status; 17% always Arrange follow-up). Three 5As factors were identified- 'Helping', 'Assessing quitting' and 'Assessing dependence'. Responses to barrier/enabler items showed greater knowledge, skills, intentions, and confidence with Assessment than Assisting; endorsement for SCS as a priority and part of midwives' professional role; and gaps in training and organisational support for SCS. Nine barrier/enabler factors were identified. Of these, the factors of 'Capability' (knowledge, skills, confidence); 'Work Environment' (service has resources, capacity, champions and values SCS) and 'Personal priority' (part of role and a priority) predicted 'Helping'. CONCLUSION: The TDF enabled systematic identification of barriers to providing SCS, and the multivariate models identified key contributors to poor implementation. Combined with qualitative data, these results have been mapped to intervention components to develop a comprehensive intervention to improve SCS.

Validity of a tool designed to assess the preventability of potentially preventable hospitalizations for chronic conditions.

BACKGROUND: Potentially preventable hospitalizations (PPH) are defined as unplanned hospital admissions which could potentially have been prevented with the provision of effective, timely outpatient care. To better understand and ultimately reduce rates of PPH, a means of identifying those which are actually preventable is required. The Preventability Assessment Tool (PAT) was designed for use by hospital clinicians to assess the preventability of unplanned admissions for chronic conditions. OBJECTIVE: The present study examined the ability of the PAT to distinguish between those unplanned admissions which are preventable and those which are not, compared to the assessments of an Expert Panel. METHODS: Data were collected between November 2014 and June 2017 at three hospitals in NSW, Australia. Participants were community-dwelling patients with unplanned hospital admissions for congestive heart failure, chronic obstructive pulmonary disease, diabetes complications or angina pectoris. A nurse and a doctor caring for the patient made assessments of the preventability of the admission using the PAT. Expert Panels made assessments of the preventability of each admission based on a comprehensive case report and consensus process. RESULTS: There was little concordance between the hospital doctors and nurses regarding the preventability of admissions, nor between the assessments of the Expert Panel and the hospital nurse or the Expert Panel and the hospital doctor. CONCLUSIONS: The PAT demonstrated poor concurrent validity and is not a valid tool for assessing the preventability of unplanned hospital admissions. The use of Expert Panels provides a more rigorous approach to assessing the preventability of such admissions.

Orofacial overgrowth with peripheral nerve enlargement and perineuriomatous pseudo-onion bulb proliferations is part of the PIK3CA-related overgrowth spectrum.

Individuals with orofacial asymmetry due to mucosal overgrowths, ipsilateral bone and dental aberrations with perineurial hyperplasia and/or perineuriomatous pseudo-onion bulb proliferations, comprise a recognizable clinical entity. In this article, we describe three individuals with this clinical entity and mosaic PIK3CA variants c.3140A>G (p. His1047Arg), c.328_330delGAA (p. Glu110del), and c.1353_1364del (p.Glu453_Leu456del). We conclude that the identification of these mosaic variants in individuals with orofacial asymmetry presenting histopathologically perineurial hyperplasia and/or intraneural pseudo-onion bulb perineurial cell proliferations supports the inclusion of this clinical entity in the PIK3CA-related overgrowth spectrum.

NAA10 polyadenylation signal variants cause syndromic microphthalmia.

BACKGROUND: A single variant in NAA10 (c.471+2T>A), the gene encoding N-acetyltransferase 10, has been associated with Lenz microphthalmia syndrome. In this study, we aimed to identify causative variants in families with syndromic X-linked microphthalmia. METHODS: Three families, including 15 affected individuals with syndromic X-linked microphthalmia, underwent analyses including linkage analysis, exome sequencing and targeted gene sequencing. The consequences of two identified variants in NAA10 were evaluated using quantitative PCR and RNAseq. RESULTS: Genetic linkage analysis in family 1 supported a candidate region on Xq27-q28, which included NAA10. Exome sequencing identified a hemizygous NAA10 polyadenylation signal (PAS) variant, chrX:153,195,397T>C, c.*43A>G, which segregated with the disease. Targeted sequencing of affected males from families 2 and 3 identified distinct NAA10 PAS variants, chrX:g.153,195,401T>C, c.*39A>G and chrX:g.153,195,400T>C, c.*40A>G. All three variants were absent from gnomAD. Quantitative PCR and RNAseq showed reduced NAA10 mRNA levels and abnormal 3' UTRs in affected individuals. Targeted sequencing of NAA10 in 376 additional affected individuals failed to identify variants in the PAS. CONCLUSION: These data show that PAS variants are the most common variant type in NAA10-associated syndromic microphthalmia, suggesting reduced RNA is the molecular mechanism by which these alterations cause microphthalmia/anophthalmia. We reviewed recognised variants in PAS associated with Mendelian disorders and identified only 23 others, indicating that NAA10 harbours more than 10% of all known PAS variants. We hypothesise that PAS in other genes harbour unrecognised pathogenic variants associated with Mendelian disorders. The systematic interrogation of PAS could improve genetic testing yields.

Missense Variants in the Histone Acetyltransferase Complex Component Gene TRRAP Cause Autism and Syndromic Intellectual Disability.

Acetylation of the lysine residues in histones and other DNA-binding proteins plays a major role in regulation of eukaryotic gene expression. This process is controlled by histone acetyltransferases (HATs/KATs) found in multiprotein complexes that are recruited to chromatin by the scaffolding subunit transformation/transcription domain-associated protein (TRRAP). TRRAP is evolutionarily conserved and is among the top five genes intolerant to missense variation. Through an international collaboration, 17 distinct de novo or apparently de novo variants were identified in TRRAP in 24 individuals. A strong genotype-phenotype correlation was observed with two distinct clinical spectra. The first is a complex, multi-systemic syndrome associated with various malformations of the brain, heart, kidneys, and genitourinary system and characterized by a wide range of intellectual functioning; a number of affected individuals have intellectual disability (ID) and markedly impaired basic life functions. Individuals with this phenotype had missense variants clustering around the c.3127G>A p.(Ala1043Thr) variant identified in five individuals. The second spectrum manifested with autism spectrum disorder (ASD) and/or ID and epilepsy. Facial dysmorphism was seen in both groups and included upslanted palpebral fissures, epicanthus, telecanthus, a wide nasal bridge and ridge, a broad and smooth philtrum, and a thin upper lip. RNA sequencing analysis of skin fibroblasts derived from affected individuals skin fibroblasts showed significant changes in the expression of several genes implicated in neuronal function and ion transport. Thus, we describe here the clinical spectrum associated with TRRAP pathogenic missense variants, and we suggest a genotype-phenotype correlation useful for clinical evaluation of the pathogenicity of the variants.

Correction: Secondary findings from clinical genomic sequencing: prevalence, patient perspectives, family history assessment, and health-care costs from a multisite study.

The originally published version of this Article contained errors in Fig. 2. The numbers below the black arrowheads were incorrect; please see incorrect Figure in associated Correction. These errors have now been corrected in the PDF and HTML versions of the Article.